arXiv:l 602.02 120v4 [cs.DS] 12 Nov 2016 


Parallel Ordered Sets Using Join 


Guy Blelloch Daniel Ferizovic 

Carnegie Mellon University Karlsruhe Institute of Technology 

guyb@cs.cmu.edu dani93.f@gmail.com 

Yihan Sun 

Carnegie Mellon University 
yihans@cs.cmu.edu 

November 15, 2016 


Abstract 

Ordered sets (and maps when data is associated with each key) are one of the most important and useful data types. 
The set-set functions union, intersection and difference are particularly useful in certain applications. Brown and 
Tarjan first described an algorithm for these functions, based on 2-3 trees, that meet the optimal O (m log ( — + 1 ) ) 
time bounds in the comparison model ( n and m < n are the input sizes). Later Adams showed very elegant algorithms 
for the functions, and others, based on weight-balanced trees. They only require a single function that is specific to the 
balancing scheme — a function that joins two balanced trees — and hence can be applied to other balancing schemes. 
Furthermore the algorithms are naturally parallel. However, in the twenty-four years since, no one has shown that the 
algorithms are work efficient (or optimal), sequential or parallel, and even for the original weight-balanced trees. 

In this paper we show that Adams’ algorithms are both work efficient and highly parallel (polylog depth) across 
four different balancing schemes — AVL trees, red-black trees, weight balanced trees and treaps. To do this we need 
careful, but simple, algorithms for JOIN that maintain certain invariants, and our proof is (mostly) generic across the 
schemes. 

To understand how the algorithms perform in practice we have also implemented them (all code except JOIN is 
generic across the balancing schemes). Interestingly the implementations on all four balancing schemes and three set 
functions perform similarly in time and speedup (more than 45x on 64 cores). We also compare the performance of 
our implementation to other existing libraries and algorithms including the standard template library (STL) imple- 
mentation of red-black trees, the multicore standard template library (MCSTL), and a recent parallel implementation 
based on weight-balanced trees. Our implementations are not as fast as the best of these on fully overlapping keys 
(but comparable), but better than all on keys with a skewed overlap (two Gaussians with different means). 
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1 Introduction 


Ordered sets and ordered maps (sets with data associated with each key) are two of the most important data types 
used in modern programming. Most programming languages either have them built in as basic types (e.g. python) or 
supply them as standard libraries (C++, C# Java, Scala, Haskell, ML). These implementations are based on some form 
of balanced tree (or tree-like) data structure and, at minimum, support lookup, insertion, and deletion in logarithmic 
time. Most also support set-set functions such as union, intersection, and difference. These functions are particularly 
useful when using parallel machines since they can support parallel bulk updates. In this paper we are interested in 
simple and efficient parallel algorithms for such set-set functions. 

The lower bound for comparison-based algorithms for union, intersection and difference for inputs of size n and 
to < n and output an ordered structur^jis log 2 ( m +") = 0(TOlog(^ + l))- Brown and Tarjan first described a 
sequential algorithm for merging that asymptotically match these bounds ED. It can be adapted for union, intersection 
and difference with the same bounds. The bound is interesting since it shows that implementing insertion with union, 
or deletion with difference, is asymptotically efficient (O(logn) time), as is taking the union of two equal sized sets 
(O(n) time). However, the Brown and Tarjan algorithm is complicated, and completely sequential. 


insert (T,fc) = 

( T L ,m,T R ) = split (T,k); 
join (T L ,k,T R ) 

delete (T, fc) = 

{TL,m,T R ) = split(T, fc); 
join2(7T , Tr) 


union(Ti,T 2 ) = 

if Ti = Leaf then T 2 

else if T 2 = Leaf then T \ 

else ( L 2 ,k 2 ,R 2 ) = expose(T 2 ); 

(L lr b,Ri) = split(Ti, fc 2 ); 

Tl = union(Li , L 2 ) || T R = union(i?i , R 2 ) ; 
join (T L ,k 2r T R ) 


split (T,k) = 

if T = Leaf then (Leaf , false. Leaf) 
else ( L,m.,R ) = expose(T); 
if k m then (L,true,i?) 
else if k < m then 

(L L ,b,L R ) = split (L,fc); 

{L L ,b, join (L R ,m,R)) 
else (R R ,b,R R ) = split (R,k); 
(join (L,m,R L ),b, R r ) 

splitLast(T) = 

( L,k,R ) = expose(T); 
if R = Leaf then (T, k) 
else (T',fc') = splitLast(R); 

( join(L, fc, T') , k') 

join2 (T l ,T r ) = 

if T r = Leaf then T R 
else ( T' L ,k ) = splitLast(Ti); 
join(T^, k,T R ) 


intersect (Ti , T 2 ) = 

if Ti = Leaf then Leaf 
else if T 2 = Leaf then Leaf 
else ( L2,k 2 ,R 2 ) = expose(T 2 ); 

(L\,b,R\) = split(Ti, fc 2 ); 

T r = intersect(Ti , L 2 ) || T R = intersect(f?i , R 2 )t 
if b = true then join(Tj,, fc 2 , T R ) 
else join2 (T L ,T R ) 

dif ference(Ti , T 2 ) = 

if Ti — Leaf then Leaf 

else if T 2 = Leaf then T 1 

else (L,2,k2 r R 2 ) = expose(T 2 ); 

(L 1 ,b,R 1 ) = split(Ti, fc 2 ); 

Tl = dif ference(Ti , L 2 ) || T R = dif f erence(i?i , R 2 ) ; 
join2 (T Lr T R ) 


Figure 1: Implementing Union, Intersect, Difference, Insert, Delete, Split, and Join2 with just Join. 
Expose returns the left tree, key, and right tree of a node. The 1 1 notation indicates the recursive calls can run in 
parallel. These are slight variants of the algorithms described by Adams m, although he did not consider parallelism. 

Adams later described very elegant algorithms for union, intersection, and difference, as well as other functions 
using weight-balanced trees, based on a single function. Join IH Hi (see Figure [l]). The algorithms are naturally 
parallel. The Join(L, fc, R) function takes a key fc and two ordered sets L and R such that L < k < R and returns the 
union of the keys ll28l|261 . Join can be used to implement Join2(L, R), which does not take the key in the middle, 
and Split(T, fc), which splits a tree at a key fc returning the two pieces and a flag indicating if fc is in T (See Section 
[4). With these three functions, union, intersection, and difference (as well as insertion, deletion and other functions) 

1 By “ordered structure” we mean any data structure that can output elements in sorted order without any comparisons. 
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are almost trivial. Because of this at least three libraries use Adams’ algorithms for their implementation of ordered 
sets and tables (Haskell 1(191 and MIT/GNU Scheme, and SML). 

JOIN can be implemented on a variety of different balanced tree schemes. Sleator and Tarjan describe an algorithm 
for Join based on splay trees which runs in amortized logarithmic time |[26l . Tarjan describes a version for red-black 
tree that runs in worst case logarithmic time |[28l . Adams describes version based on weight-balanced trees mQ 
Adams’ algorithms were proposed in an international competition for the Standard ML community, which is about 
implementations on “set of integers”. Prizes were awarded in two categories: fastest algorithm, and most elegant yet 
still efficient program. Adams won the elegance award, while his algorithm is as fast as the fastest program for very 
large sets, and was faster for smaller sets. Adams’ algorithms actually show that in principle all balance criteria for 
search trees can be captured by a single function JOIN, although he only considered weight-balanced trees. 

Surprisingly, however, there have been almost no results on bounding the work (time) of Adams’ algorithms, in 
general nor on specific trees. Adams informally argues that his algorithms take 0(n + to) work for weight-balanced 
tree, but that is a very loose bound. Blelloch and Miller later show that similar algorithms for treaps (6), are optimal for 
work (i.e. 0 (to log ( — + l) ) ), and also parallel. Their algorithms, however, are specific for treaps. The problem with 
bounding the work of Adams’ algorithms, is that just bounding the time of Split, Join and Join2 with logarithmic 
costs is not sufficient^] One needs additional properties of the trees. 

The contribution of this paper is to give first work-optimal bounds for Adams’ algorithms. We do this not only for 
the weight-balanced trees, we bound the work and depth of Adams’ algorithms (union, intersection and difference) 
for four different balancing shemes: AVL trees, red-black trees, weight-balanced trees and treaps. We analyze exactly 
the algorithms in Figure [I] and the bounds hold when either input tree is larger. We show that with appropriate (and 
simple) implementations of JOIN for each of the four balancing schemes, we achieve asymptotically optimal bounds 
on work. Furthermore the algorithms have O(lognlogm) depth, and hence are highly parallel. To prove the bounds 
on work we show that our implementations of JOIN satisfy certain conditions based on a rank we define for each tree 
type. In particular the cost of JOIN must be proportional to the difference in ranks of two trees, and the rank of the 
result of a join must be at most one more than the maximum rank of the two arguments. 

In addition to the theoretical analysis of the algorithms, we implemented parallel versions of all of the algorithms 
on all four tree types and describe a set of experiments. Our implementation is generic in the sense that we use 
common code for the algorithms in Figure [I] and only wrote specialized code for each tree type for the JOIN function. 
Our implementations of JOIN are as described in this paper. We compare performance across a variety of parameters. 
We compare across the tree types, and interestingly all four balance criteria have very similar performance. We 
measure the speedup on up to 64 cores and achieve close to a 46-fold speedup. We compare to other implementations, 
including the set implementation in the C++ Standard Template library (STL) for sequential performance, and parallel 
weight-balanced B-trees (WBB-trees) fT2l and the multi-core standard template library (MCSTL) lfl3l for parallel 
performance. We also compare for different data distributions. The conclusion from the experiments is that although 
not always as fast as (WBB-trees) m on uniform distributions, the generic code is quite competitive, and on keys 
with a skewed overlap (two Gaussians with different means), our implementation is much better than all the other 
baselines. 

Related Work Parallel set operations on two ordered sets have been well-studied, but each previous algorithm only 
works on one type of balance tree. Paul, Vishkin, and Wagener studied bulk insertion and deletion on 2-3 trees in the 
PRAM model ll24l . Park and Park showed similar results for red-black trees |23j . These algorithms are not based 
on Join. Katajainen lfl~6ll claimed an algorithm with 0(?7ilog(^ + 1)) work and O(logn) depth using 2-3 tree, but 
was shown to contain some bugs so the bounds do not hold |6). Akhremtsev and Sanders j4j (unpublished) recently 
fixed that and proposed an algorithm based on (a, 6)-trees with optimal work and 0(log n) depth. Their algorithm 
only works for (a, /i) -trees, and they only gave algorithms on UNION. Besides, our algorithm is much simpler and 
easy to be implemented than theirs. Blelloch and Miller showed a similar algorithm as Adams’ (as well as ours) 
on treaps with optimal work and O(logn) depth on a EREW PRAM with scan operation using pipelines (implying 
0(log n log to) depth on a plain EREW PRAM, and 0(log* m log n) depth on a plain CRCW PRAM). The pipelining 

2 Adams’ version had some bugs in maintaining the balance, but these were later fixed (1411271. 

^Bounding the cost of JOIN, SPLIT, and JOIN2 by the logarithm of the smaller tree is probably sufficient, but implementing a data structure with 
such bounds is very much more complicated. 
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is quite complicated. Our focus in this paper is in showing that very simple algorithm are work efficient and have 
polylogarithmic depth, and less with optimizing the depth. 

Many researchers have considered concurrent implementations of balanced search trees (e.g., 02 HSUS ED). 
None of these is work efficient for Union since it is necessary to insert one tree into the other taking work at least 
0(m log n). Furthermore the overhead of concurrency is likely to be very high. 

2 Preliminaries 

A binary tree is either a Leaf, or a node consisting of a left binary tree Tr, a value (or key) v, and a right binary 
tree Tr, and denoted Node(Tr,, v, Tr). The size of a binary tree, or |T|, is 0 for a Leaf and \Tr\ + \Tr\ + 1 for a 
Node(Ti, v, Tr). The weight of a binary tree, or w(T ), is one more than its size (i.e., the number of leaves in the 
tree). The height of a binary tree, or h(T), is 0 for a Leaf, and max(/i(T/,), Ii(Tr)) + 1 for a Node(Tf,, v,Tr). 
Parent , child, ancestor and descendant are defined as usual (ancestor and descendant are inclusive of the node itself). 

The left spine of a binary tree is the path of nodes from the root to a leaf always following the left tree, and the right 
spine the path to a leaf following the right tree. The in-order values of a binary tree is the sequence of values returned 
by an in-order traversal of the tree. 

A balancing scheme for binary trees is an invariant (or set of invariants) that is true for every node of a tree, and 
is for the purpose of keeping the tree nearly balanced. In this paper we consider four balancing schemes that ensure 
the height of every tree of size n is bounded by Oflog n). For each balancing scheme we define the rank of a tree, or 

r(T). 

AVL trees 0 have the invariant that for every Node (7 7,. v, Tr), the height of Tr and Tr differ by at most one. 

This property implies that any AVL tree of size n has height at most log^f n + 1), where f is the golden ratio. 

For AVL trees r(T) = h(T) - 1. 

Red-black (RB) trees 0 associate a color with every node and maintain two invariants: (the red rule) no red 
node has a red child, and (the black rule) the number of black nodes on every path from the root down to a leaf is 
equal. Unlike some other presentations, we do not require that the root of a tree is black. Our proof of the work 
bounds requires allowing a red root. We define the black height of a node T, denoted h(T) to be the number of black 
nodes on a downward path from the node to a leaf (inclusive of the node). Any RB tree of size n has height at most 
2 log 2 (n + 1). In RB trees r{T) = 2 (h{T) - 1) if T is black and r(T) = 2 h(T) - 1 if T is red. 

Weight-balanced (WB) trees with parameter a (also called BB[a] trees) ll22l maintain for every T = Node(7V, , v, Tr) 
the invariant a < — 1 — a - We say two weight-balanced trees Ti and T% have like weights if Node(7i, v, Tf) is 

weight balanced. Any a weight-balanced tree of size n has height at most log ^ 1 n. For < a < 1 — insertion 
and deletion can be implemented on weight balanced trees using just single and double rotations Il22l l7ll. We require 
the same condition for our implementation of Join, and in particular use a = 0.29 in experiments. For WB trees 
r(T) = riog 2 (w;(T))l - 1 . 

Treaps Ii25l associate a uniformly random priority with every node and maintain the invariant that the priority 
at each node is no greater than the priority of its two children. Any treap of size n has height Oflog n) with high 
probability (w.h.pjj For treaps r(T) = [log 2 (w(T))] — 1. 

For all the four balancing schemes r(T) = 0(log(|Tj + 1)). The notation we use for binary trees is summarized 
in Table Q] 

A Binary Search Tree (BST) is a binary tree in which each value is a key taken from a total order, and for which 
the in-order values are sorted. A balanced BST is a BST maintained with a balancing scheme, and is an efficient way 
to represent ordered sets. 

Our algorithms are based on nested parallelism with nested fork-join constructs and no other synchronization or 
communication among parallel tasks All algorithms are deterministic. We use work (W) and span (S) to analyze 
asymptotic costs, where the work is the total number of operations and span is the critical path. We use the simple 
composition rules W(e 1 || ef) = W(ei) + W(e 2 ) + 1 and S(e 1 || 62 ) = max(S(ei), Sfa)) + 1. For sequential 

4 Here w.h.p. means that height Ofelog n) with probability at least ! — 1 /n r (c is a constant) 

' This does not preclude using our algorithms in a concurrent setting. 


4 


Notation 

Description 

|T| 

The size of tree T 

h{T) 

The height of tree T 

h(T) 

The black height of an RB tree T 

r{T) 

The rank of tree T 

w{T) 

The weight of tree T (i.e, |T| + 1) 

P(T) 

The parent of node T 

k(T) 

The value (or key) of node T 

L{T) 

The left child of node T 

R{T) 

The right child of node T 

expos e(T) 

c L(T),k(T),R(T )) 


Table 1 : Summary of notation. 


computation both work and span compose with addition. Any computation with W work and S span will run in time 
T < + S assuming a PRAM (random access shared memory) with P processors and a greedy scheduler JU®. 


3 The JOIN Function 

Here we describe algorithms for JOIN for the four balancing schemes we defined in Section[2] The function JOIN (Tr, k. Tr) 
takes two binary trees Tr and Tr, and a value k, and returns a new binary tree for which the in-order values are a 
concatenation of the in-order values of Tr, then k, and then the in-order values of Tr. 

As mentioned in the introduction and shown in Section |4j JOIN fully captures what is required to rebalance a tree 
and can be used as the only function that knows about and maintains the balance invariants. For AVL, RB and WB 
trees we show that JOIN takes work that is proportional to the difference in rank of the two trees. For treaps the work 
depends on the priority of k. All versions of JOIN are sequential so the span is equal to the work. Due to space 
limitations, we describe the algorithms, state the theorems for all balancing schemes, but only show a proof outline for 
AVL trees. 


1 joinRight (T L ,k,T R ) = 

2 ( l,k',c ) = expose(Ty); 

3 if h(c) < h(T R ) + 1 then 

4 T' = Node(c, k, Tr) ; 

5 if h(T') < h(l) + 1 then Node(l, fc\ T’) 

6 else rotateLeft(Node(l, fc', rotateRight(T'))) 

7 else 

8 T 1 = joinRight (c, k, Tr) ; 

9 T" = Nod e(l,k’,T'); 

10 if h{T') < Ml) + 1 then T" 

11 else rotateLeft(T") 

12 join(T £/ ,fc,T R ) = 

13 if h(T R ) > h(T R ) + 1 then joinRight (TR,k,T R ) 

14 else if h(TR) > h(T R ) + 1 then joinLeft(Tj,, fc, Tr) 

15 else Node(2T, k, Tr) 


Figure 2: AVL JOIN algorithm. 

AVL trees. Pseudocode for AVL JOIN is given in Figure [2] and illustrated in Figure [6] Every node stores its own 
height so that h(-) takes constant time. If the two trees Tr and Tr differ by height at most one, JOIN can simply 
create a new Nocie(7V, Tr). However if they differ by more than one then rebalancing is required. Suppose that 
h(TR) > /i(Tr) + 1 (the other case is symmetric). The idea is to follow the right spine of Tr until a node c for which 
h(c) < h{TR) + 1 is found (line [3|l. At this point a new Node(c, k, Tr) is created to replace c (line [4]). Since either 
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1 joinRightRB(Tx / , k,T R ) = 

2 if (r{T L ) = [r(Tfl)/2j x 2) then 

3 Node(Ti, (k, red), Tr); 

4 else 

5 (L', (k' , c'), -R')=expose(T£ / ); 

6 T' = Node(Z/ / , (A:', c 7 ) , joinRightRB(fJ', fc, T/j)); 

7 if (c'=black) and (c(f?(T')) = c(R(R(T')))=red) then 

8 c(R(i?(T')))=black; 

9 T"=rotateLeft(T') 

10 else T" 

11 joinRB(TR , fc, Tr) = 

12 if |r(T L )/2j > |t(Tr)/ 2J then 

13 T' = joinRightRB(TR, fc, Tr); 

14 if (c(T')=red) and (c(-R(T , ))=red) then 

15 Node(L(T , ),(fc(T'),black),.R(T')) 

16 else T' 

17 else if \r(T R )/2\ > \r(T L )/2\ then 

18 T' = joinLeftRB(T£,, k, Tr); 

19 if (c(T')=red) and (c(L(T')) = red) then 

20 Node(L(T'), (k(T'), black.) , R{T')) 

21 else T' 

22 else if (c(7Y)=black) and (c(T/j)=black) then 

23 Node(Ti, (fc, red), Tr) 

24 else NodefTj,, (k, black), Tr) 


Figure 3: RB JOIN algorithm. 


1 joinRightWB(TR, k,T R ) = 

2 (l, k' , c)=expose(Ti,) ; 

3 if (balance(|Tx,|, |Tr|) then Node(Tx,, k, Tr)); 

4 else 

5 T' = joinRightWB(c, Ai, Tr); 

6 = expose(T'); 

7 if like(|/|, |T'|) then Nod e(l,k',T') 

8 else if (like(|Z|, |Zi|)) and (like(|Z| + \l± |, ri)) then 

9 rotateLeft(Node(Z, k' , T')) 

10 else rotateLeft(Node(Z, k' , rotateRight(T / ))) 

11 joinWB(Tx,, A:, Tr) = 

12 if heavy(T£ / , Tr) then joinRightWB(TR, k, Tr) 

13 else if heavy (T r,Tr) then joinLeftWB(Tx / , A;, Tr) 

14 else Node(T£, k, Tr) 


Figure 4: WB Join algorithm. 


1 joinTreap(Tx / , k, Tr) = 

2 if prior(A), Aii) and prior(Ai, A 12 ) then Node(TR, k, Tr) 

3 else (Zi, fci,ri)=expose(T£,); 

4 (Z 2 , &2, r2)=expose(T^); 

5 if prior(fci, A 12 ) then 

6 Node(Zi , Aii r joinTreap(r*i, k, Tr)) 

7 else Node( joinTreap(Tt, fc, I 2 ), L 2 , r^) 


Figure 5: Treap JOIN algorithm. 


h(c) = h(Tfi) or h(c) = h{T R ) + 1, the new node satisfies the AVL invariant, and its height is one greater than c. The 
increase in height can increase the height of its ancestors, possibly invalidating the AVL invariant of those nodes. This 
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can be fixed either with a double rotation if invalid at the parent (line|6| or a single left rotation if invalid higher in the 
tree (linefTTj), in both cases restoring the height for any further ancestor nodes. The algorithm will therefore require at 
most two rotations. 



Figure 6: An example for Join on AVL trees (/i(Tr) > h(Tn) + 1). We first follow the right spine of Tr until a 
subtree of height at most h{T r ) + 1 is found (i.e., T-± rooted at c). Then a new Node(c, Ic,Tr) is created, replacing 
c (Step 1). If h(T-\ ) = h and h(Tf) = h + 1, the node p will no longer satisfy the AVL invariant. A double rotation 
(Step 2) restores both balance and its original height. 


Lemma 1. For two AVL trees Tr and Tr, the AVL JOIN algorithm works correctly, runs with 0(|/i(Tr) — 
work, and returns a tree satisfying the AVL invariant with height at most 1 + max(/i(Ti), h(Tp)). 

Proof outline. Since the algorithm only visits nodes on the path from the root to c, and only requires at most two 
rotations, it does work proportional to the path length. The path length is no more than the difference in height of the 
two trees since the height of each consecutive node along the right spine of Tr differs by at least one. Along with the 
case when H{Tr) > /i(Tr) + 1, which is symmetric, this gives the stated work bounds. The resulting tree satisfies the 
AVL invariants since rotations are used to restore the invariant (details left out). The height of any node can increase 
by at most one, so the height of the whole tree can increase by at most one. □ 

Red-black Trees. Tarjan describes how to implement the JOIN function for red-black trees l28l . Here we describe 
a variant that does not assume the roots are black (this is to bound the increase in rank by Union). The pseudocode 
is given in Figure [ 3 ] We store at every node its black height /<(■). The first case is when Ii(Tr) = h(Tjf). Then if 
both ^(Tr) and k{T£) are black, we create red Node(TR, k , Tr), otherwise we create black Node(TR, fc, Tr). The 
second case is when /i(Tr) < /i(Tr) = h (the third case is symmetric). Similarly to AVL trees. Join follows the right 
spine of Tr until it finds a black node c for which h{c) = /i(Tr). It then creates a new red Node(c, k , Tr) to replace 
c. Since both c and Tr have the same height, the only invariant that can be violated is the red rule on the root of Tr, 
the new node, and its parent, which can all be red. In the worst case we may have three red nodes in a row. This is 
fixed by a single left rotation: if a black node v has R{v ) and R(R(v)) both red, we turn R(R(v)) black and perform 
a single left rotation on v. The update is illustrated in Figure [7] The rotation, however can again violate the red rule 
between the root of the rotated tree and its parent, requiring another rotation. The double-red issue might proceed up 
to the root of Tr. If the original root of Tr is red, the algorithm may end up with a red root with a red child, in which 
case the root will be turned black, turning Tr rank from 2/t 1 to 2 h. If the original root of Tr is black, the algorithm 

may end up with a red root with two black children, turning the rank of Tr from 2h — 2 to 2h — 1. In both cases the 
rank of the result tree is at most 1 + 7'(Tr). 

Lemma 2. For two RB trees Tr and Tr, the RB JOIN algorithm works correctly, runs with 0(\r(Tif) — r(TR)|) work, 
and returns a tree satisfying the red-black invariants and with rank at most 1 + max(r(TR), t(Tr)). 

The proof is similar as Lemma[T| 

Weight Balanced Trees. We store the weight of each subtree at every node. The algorithm for joining two weight- 
balanced trees is similar to that of AVL trees and RB trees. The pseudocode is shown in Figure[4] The like function 
in the code returns true if the two input tree sizes are balanced, and false otherwise. If Tr and Tr have like weights 
the algorithm returns a new Node(TR, k , Tr). Suppose |Tr| < |Tr|, the algorithm follows the right branch of Tr 
until it reaches a node c with like weight to Tr. It then creates a new Node(c, k. Tr) replacing c. The new node will 
have weight greater than c and therefore could imbalance the weight of c’s ancestors. This can be fixed with a single 
or double rotation (as shown in Figure [8} at each node assuming a is within the bounds given in Section[2] 
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Figure 7: An example of JOIN on red-black trees ( h = h( r Frk) > Ii(Tr)). We follow the right spine of Tr until we 
find a black node with the same black height as Tr (i.e., c). Then a new red Node(c, k. Tr) is created, replacing c 
(Step 1). The only invariant that can be violated is when either c’s previous parent p or Tr’s root d is red. If so, a left 
rotation is performed at some black node. Step 2 shows the rebalance when p is red. The black height of the rotated 
subtree (now rooted at p) is the same as before ( h + 1), but the parent of p might be red, requiring another rotation. If 
the red-rule violation propagates to the root, the root is either colored red, or rotated left (Step 3). 

Lemma 3. For two a weight-balanced trees Tl and Tr and a < 1 — ^ ss 0.29, the weight-balanced JOIN algorithm 
works correctly, runs with 0( \ log(w(Ti,)/w(T}i))\) work, and returns a tree satisfying the a weight-balance invariant 
and with rank at most 1 + ma x(r(T£,), t'(T'r)). 

The proof is shown in the Appendix. 

Treaps. The treap Join algorithm (as in Figure [5j first picks the key with the highest priority among k, k{Ti) and 
k(Tjf) as the root. If k is the root then the we can return Node(T7,, k, Tr). Otherwise, WLOG, assume fc(Tj,) has 
a higher priority. In this case A: (7),) will be the root of the result, A (77, ) will be the left tree, and R(Tr), k and Tr 
will form the right tree. Thus Join recursively calls itself on R(Tr), k and Tr and uses result as fc(77,)’s right child. 
When k(TR) has a higher priority the case is symmetric. The cost of JOIN is therefore the depth of the key k in the 
resulting tree (each recursive call pushes it down one level). In treaps the shape of the result tree, and hence the depth 
of k, depend only on the keys and priorities and not the history. Specifically, if a key has the I th highest priority among 
the keys, then its expected depth in a treap is O(logt) (also w.h.p.). If it is the highest priority, for example, then it 
remains at the root. 

Lemma 4 . For two treaps Tr and Tr, if the priority ofk is the t-th highest among all keys in Tr U {k} U Tr, the treap 
JOIN algorithm works correctly, runs with 0(\ogt + 1) work in expectation and w.h.p., and returns a tree satisfying 
the treap invariant with rank at most 1 + max(r(T^), r(Tfj)). 

From the above lemmas we can get the following fact for Join. 

Theorem 1. For AVL, RB and WB trees }oin(Tr, k, Tr) does 0(|r(T7,) — r(T7j)|) work. For treaps Join does 
0(log t) work in expectation ifk has the t-th highest priority among all keys. For AVL, RB, WB trees and treaps, JOIN 
returns a tree T for which the rank satisfies ma x(r(X7,), t{Tr)) < r(T) < 1 + max(r(T£,), t(Tr)). 
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Figure 8: An illustration of single and double rotations possibly needed to rebalance weight-balanced trees. In the 
figure the subtree rooted at u has become heavier due to joining in Tr and its parent v now violates the balance 
invariant. 


4 Other Functions Using JOIN 

In this section, we describe algorithms for various functions that use just Join. The algorithms are generic across 
balancing schemes. The pseudocodes for the algorithms in this section is shown in Figure|T] 

Split. For a BST T and key k, Split(T, k) returns a triple (7/,. b, Tr), where Tr (Tr) is a tree containing all keys 
in T that are less (larger) than k, and b is a flag indicating whether k £ T. The algorithm first searches for k in T, 
splitting the tree along the path into three parts: keys to the left of the path, k itself (if it exists), and keys to the right. 
Then by applying JOIN, the algorithm merges all the subtrees on the left side (using keys on the path as intermediate 
nodes) from bottom to top to form Tr, and merges the right parts to form Tr. Figure [9] gives an example. 


Split T with key 42: 




Figure 9: An example of Split in a BST with key 42. We first search for 42 in the tree and split the tree by the 
searching path, then use Join to combine trees on the left and on the right respectively, bottom-top. 

Theorem 2. The work of Split (T, k) is 0(log|T|) for all balancing schemes described in Section [i] (w.h.p. for 
treaps). The two resulting trees Tr and Tr will have rank at most r(T). 
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Proof. We only consider the work of joining all subtrees on the left side. The other side is symmetric. Suppose we have 
l subtrees on the left side, denoted from bottom to top as Xj, T 2 , . . . Tj. We have that r(Xj) < r(T 2 ) < • • • < r(T{). 
As stated above, we consecutively join Tj and T 2 returning T 2 , then join 7’j with 7 3 returning 7'j and so forth, until 
all trees are merged. The overall work of Split is the sum of the cost of l — 1 JOIN functions. 

For AVL trees, red-black trees and weight-balanced trees, recall Theorem [I] that we have r(T') < r{Tf) + 1, so 
r(T') < r(Tf) + 1 < r(7j + 1 ) + 1. According to Theorem[l] the work of a single operation is 0(|r(7j + i )-r(T>)\). 
The overall complexity is \ r (R+i) - r ( T i)\ < E!:=i r ( T i+ 1 ) “ r ( T l) + 2 = 0(r(T )) = 0(log |T|). 

For treaps, each JOIN uses the key with the highest priority since the key is always on a upper level. Hence by 
Lemma[4] the complexity of each JOIN is 0(1) and the work of split is at most 0(log |T|) w.h.p. 

Also notice that when getting the final result 7’/-, and Tr, the last step is a Join on two trees, the larger one of 
which is a subtree of the original T. Thus the rank of the two trees to be joined is of rank at most r(T) — 1, according 
to Theorem[l]we have t(Tr) and t{Tr) at most r(T). □ 

Join2. Join2(7’/\. Tr) returns a binary tree for which the in-order values is the concatenation of the in-order values 
of the binary trees 7) and Tr (the same as Join but without the middle key). For BSTs, all keys in 7) have to be less 
than keys in Tr. Join 2 first finds the last element k (by following the right spine) in If and on the way back to root, 
joins the subtrees along the path, which is similar to Split Tr by k. We denote the result of dropping k in Tr as T' L . 
Then .1 0IN(7’( , k. Tr) is the result of Join 2. Unlike Join, the work of J0IN2 is proportional to the rank of both trees 
since both Split and Join take at most logarithmic work. 

Theorem 3. The work of Join2(T£,, Tr) is 0{t{Tr) + t{Tr)) for all balancing schemes described in Section [i] 
(bounds are w.h.p for treaps). 

Union, Intersect and Difference. Union(7i, T 2 ) takes two BSTs and returns a BST that contains the union of all 
keys. The algorithm uses a classic divide-and-conquer strategy, which is parallel. At each level of recursion, 7) is split 
by k(T 2 ), breaking 7) into three parts: one with all keys smaller than k(flV) (denoted as L \ ), one in the middle either 
of only one key equal to k(T 2 ) (when k(T 2 ) £ Tf) or empty (when k(T 2 ) ^ Tf), the third one with all keys larger than 
k(T 2 ) (denoted as R\ ). Then two recursive calls to Union are made in parallel. One unions L(T 2 ) with L \ , returning 
Tr, and the other one unions R(T 2 ) with R\, returning Tr. Finally the algorithm returns Jo IN (7)., k(T 2 ),TR), which 
is valid since k(T 2 ) is greater than all keys in Tr are less than all keys in Tr. 

The functions Intersect(Xi, T 2 ) and DIFFERENCE^, T 2 ) take the intersection and difference of the keys in 
their sets, respectively. The algorithms are similar to UNION in that they use one tree to split the other. However, the 
method for joining is different. For INTERSECT, Join 2 is used instead of Join if the root of the first is not found in 
the second. For DIFFERENCE, Join 2 is used anyway because k(T 2 ) should be excluded in the result tree. The base 
cases are also different in the obvious way. 

Theorem 4 (Main Theorem). For all the four balance schemes mentioned in Section [77] the work and span of the 
algorithm (as shown in Figure^ 7]) of UNION, INTERSECT or DIFFERENCE on two balanced BSTs of sizes m and n 
(n > m) is O^mlog^ 1 - l)) and O(lognlogm) respectively (the bound is in expectation for treaps). 

A generic proof of Theorem[4]working for all the four balancing schemes will be shown in the next section. 

The work bound for these algorithms is optimal in the comparison-based model. In particular considering all possi- 
ble interleavings, the minimum number of comparisons required to distinguish them is log = © (jn log ( — + l) ) 

02 . 

Other Functions. Many other functions can be implemented with Join. Pseudocode for Insert and Delete was 
given in Figure[l] For a tree of size n they both take (7 (log n) work. 


5 The Proof of the Main Theorem 

In this section we prove Theorem[4] for all the four balance schemes (AVL trees, RB trees, WB trees and treaps) and all 
three set algorithms (UNION, INTERSECT, DIFFERENCE) from Figure[l] For this purpose we make two observations. 
The first is that all the work for the algorithms can be accounted for within a constant factor by considering just the 
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Notation 

Description 

T P 

The pivot tree 

T d 

The decomposed tree 

n 

max(|T p |, \T d \) 

m 

min(|T p |, \T d \) 

T p (v),v £ T p 

The subtree rooted at v in T p 

T d {v),v £ T p 

The tree from 7)/ that v splits]^] 

Si 

The number of nodes in layer i 

v kj 

The j-th node on layer k in T p 

d(v ) 

The number of nodes attached to a layer 
root v in a treap 


Table 2: Descriptions of notations used in the proof. 


work done by the Splits and the Joins (or Join2s), which we refer to as split work and join work , respectively. This 
is because the work done between each split and join is constant. The second observation is that the split work is 
identical among the three set algorithms. This is because the control flow of the three algorithms is the same on the 
way down the recursion when doing Splits — the algorithms only differ in what they do at the base case and on the 
way up the recursion when they join. 

Given these two observations, we prove the bounds on work by first showing that the join work is asymptotically 
at most as large as the split work (by showing that this is true at every node of the recursion for all three algorithms), 
and then showing that the split work for UNION (and hence the others) satisfies our claimed bounds. 

We start with some notation, which is summarized in Table [2] In the three algorithms the first tree ( T) ) is split 
by the keys in the second tree (T 2 ). We therefore call the first tree the decomposed tree and the second the pivot tree , 
denoted as Td and T p respectively. The tree that is returned is denoted as T r . Since our proof works for either tree 
being larger, we use m = min(|T p |, \Td\) and n = max(|T p |, \Td\). We denote the subtree rooted at v £ T p as T p (v), 
and the tree of keys from Td that v splits as Td(v) (i.e., Split(u, Td(v)) is called at some point in the algorithm). For 
v £ T p , we refer to \Td(v)\ as its splitting size. 

Figure[lO](a) illustrates the pivot tree with the splitting size annotated on each node. Since Split has logarithmic 
work, we have, 

split work = O j (log \Td(v)\ + 1) ) , 

\UGTjj J 

which we analyze in Theorem[6] We first, however, show that the join work is bounded by the split work. We use the 
following Lemma, which is proven in the appendix. 

Lemma 5. For T r =Union(T p , Td) on AVL, RB or WB trees, ifr(T p ) > r(Td) then r[T r ) < r(T p ) + r(Td). 

Theorem 5. For each function call to UNION, INTERSECT or DIFFERENCE on trees T p (v) and Td(v), the work to do 
the JOIN (or Join2) is asymptotically no more than the work to do the SPLIT. 

Proof. For Intersect or Difference, the cost of Join (or Join2) is 0(log(|r r | + 1)). Notice that Differ- 
ence returns the keys in Td\T p . Thus for both INTERSECT and DIFFERENCE we have T r C Td- The join work is 
0(log(|T r | + 1)), which is no more than 0(log(|7d| + 1)) (the split work). 

For Union, if r(T p ) < r(Td), the Join will cost 0{r{Td)), which is no more than the split work. 

Consider r(T p ) > r(7j/) for AVL, RB or WB trees. The rank of L(T p ) and R(T p ), which are used in the recursive 
calls, are at least r(T p ) — 1. Using Lemma Bj the rank of the two trees returned by the two recursive calls will be at 
least ( r(T p ) — 1) and at most ( r(T p ) + r(Tdf), and differ by at most 0(r(Td)) = 0(log \Td\ + 1). Thus the join cost 
is 0(log \Td\ + 1), which is no more than the split work. 

6 The nodes in T,/ (v) form a subset of T,/. but not necessarily a subtree. See details later. 
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Figure 10: An illustration of splitting tree and layers. The tree in (a) is T p , the dashed circle are the exterior nodes. 
The numbers on the nodes are the sizes of the tree from Td to be split by this node, i.e., the “splitting size” 7j/(+)j. In 
(b) is an illustration of layers on an AVL tree. 


Consider r{T p ) > r{Td) for treaps. If r(T p ) > r(Td), then \T p \ > |T<j|. The root of T p has the highest priority 
among all \T p \ keys, so on expectation it takes at most the < 2-th place among all the |T^| + \T p \ keys. From 

Lemma[4]we know that the cost on expectation is E[log t] + 1 < log E[f] + 1 < log 2 + 1, which is a constant. □ 

This implies the total join work is asymptotically bounded by the split work. 

We now analyze the split work. We do this by layering the pivot tree starting at the leaves and going to the root 
and such that nodes in a layer are not ancestors of each other. We define layers based on the ranks and denote the size 
of layer i as s,. We show that s, shrinks geometrically, which helps us prove our bound on the split work. For AVL 
and RB trees, we group the nodes with rank i in layer i. For WB trees and treaps, we put a node v in layer i iff v has 
rank i and v’s parent has rank strictly greater than i. Figure |T()](b) shows an example of the layers of an AVL tree. 

Definition 1. In a BST, a set of nodes V is called a disjoint set if and only if for any two nodes V\, V 2 in V , V\ is not 
the ancestor of v 2 . 

Lemma 6. For any disjoint set V C T p , Y^ v ev I Td(v)\ < \Td\. 

The proof of this Lemma is straightforward. 

Lemma 7. For an AVL, RB, WB tree or a treap of size N, each layer is a disjoint set, and Si < — holds for some 
constant c > 1. 

Proof For AVL, RB, WB trees and treaps, a layer is obviously a disjoint set: a node and its ancestor cannot lie in the 
same layer. 

For AVL trees, consider a node in layer 2, it must have at least two descendants in layer 0. Thus s 0 > 2s2. Since 
an AVL tree with its leaves removed is still an AVL tree, we have Si > 2sj + 2- Since So an d si are no more than N, 
we can get that Sj < 9[ ^ 2J . 

For RB trees, the number of black nodes in layer 2z is more than twice as many as in layer 2 {i + 1) and less than 
four times as many as in layer 2 (i + 1), i.e., S 21 > 2.s'2,;+2- Also, the number of red nodes in layer 2/ + 1 is no more 
than the black nodes in layer 2 i. Since sq and Si are no more than N, Si < 2l j^ 2J . 
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For WB trees and treaps, the rank is defined as [log 2 ('u;(7’))] — 1, which means that a node in layer i has weight 
at least 2 i + 1. Thus s* < (N + l)/(2 i + 1) < N/2\ □ 

Not all nodes in a WB tree or a treap are assigned to a layer. We call a node a layer root if it is in a layer. We attach 
each node u in the tree to the layer root that is u’s ancestor and has the same rank as u. We denote d(v) as the number 
of descendants attached to a layer root v. 

Lemma 8. For WB trees and treaps, if v is a layer root, d(v) is less than a constant (in expectation for treaps). 
Furthermore, the random variables d(v) for all layer roots in a treap are i.i.d. (See the proof in the Appendix.) 

By applying Lemma [7] and [8] we prove the split work. In the following proof, we denote Vkj as the j-th node in 
layer k. 

Theorem 6. The split work in UNION, INTERSECT and DIFFERENCE on two trees of size to and n is 0(m log(^ + l)). 

Proof. The total work of SPLIT is the sum of the log of all the splitting sizes on the pivot tree O ^^ ugT log(|Td(u)| + 1) 
Denote l as the number of layers in the tree. Also, notice that in the pivot tree, in each layer there are at most \Td\ 
nodes with \Td(vkj)\ > 0. Since those nodes with splitting sizes of 0 will not cost any work, we can assume s i < \T d \. 
We calculate the dominant term YhveT log(|2d(u)| + 1) i n the complexity by summing the work across layers: 


l Sk l 

^^iogdTrf^OI+i) < ^ Sfc log 
k = 0 j—1 k = 0 


Y.Mv k j)\ + i 


Sk 


y>iog(H+i 

^ \ s k 


k = 0 


We split it into two cases. If \T d \ > \T p \, ^ always dominates 1. we have: 

1 ( 

^Sfcbg - 


\ ^ , ( \Td\ 

L^ log — 

k = 0 V Sk 


1 = 


k=0 

l 


A~ + 1 
\Sk 


< E 


< 


< 


C lfe/2J 

fc=0 

1/2 

\ — > 771 

2 E^ lo g 

k=0 
1/2 


log 


TO / C [fc/2J 


m/c k 


1/2 


_ v — "\ TIL . 'll v — "V 771 

2 E “fc log— + 2 E fc ^ 

k = o fc=r 

0\m\og — ) + O(m) 

\ w/ 


= 0[m 


(' 


TO, 

n 


+ 1 


(1) 

( 2 ) 


(3) 


If \T d \ < \T p \. T f‘ can be less than 1 when k is smaller, thus the sum should be divided into two parts. Also note 
that we only sum over the nodes with splitting size larger than 0. Even though there could be more than \T d \ nodes in 
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one layer in T p , only \Td\ of them should count. Thus we assume Sk < |T^|, and we have: 
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From (|TJ» to <[2jl and (|4j> to (|5j> we apply Lemma[7]and the fact that fix) = x log(^ + 1) is monotonically increasing 
when x < n. 

For WB trees and treaps, the calculation above only includes the log of splitting size on layer roots. We need to 
further prove the total sum of the log of all splitting size is still O (to log ( + l) ) . Applying Lemma[8j the expectation 
is less than: 
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l Xk 

^EE d{v kj )log((T d (v kj ) + 1) 

k—0 j—1 


l Xk 


= E [d(v kj )] x 2y]y^log ({T d (v kj ) + 1) 

k—0 j—l 


°(mlo g (X + 1)) 


For WB trees d(v k j) is no more than a constant, so we can also come to the same bound. 

To conclude, the split work on all four balancing schemes of all three functions is O ( to log ( ^ + l) ) . □ 

Theorem 7. The total work of Union, INTERSECT or DIFFERENCE of all four balancing schemes on two trees of 
size to and n(m> n) is 0(mlog( ^ + l)). 

This directly follows Theorem [5] and [6] 

Theorem 8. The span of Union and INTERSECT or DIFFERENCE on all four balancing schemes is Otfogn logm). 
Here n and m are the sizes of the tw’o tree. 

Proof For the span of these algorithms, we denote D(h\ , / 12 ) as the span on UNION, INTERSECT or DIFFERENCE on 
two trees of height hi and h 2 . According to Theorem [5] the work (span) of Split and Join are both 0(log 7)/|) = 
0(h{T d )). We have: 

D(h(T p ),h(T d )) < D(h(T p ) - 1, h(T d )) + 2 h(T d ) 

Thus D(h(T p ), h{T d )) < 2 h(T p )h{T d ) = O(lognlogm). 0 

Combine Theorem [7] and [8] we come to Theorem |4] 
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6 Experiments 


To evaluate the performance of our algorithms we performed several experiments across the four balancing schemes 
using different set functions, while varying the core count and tree sizes. We also compare the performance of our 
implementation to other existing libraries and algorithms. 

Experiment setups and baseline algorithms For the experiments we use a 64-core machine with 4 x AMD Opteron(tm) 
Processor 6278 (16 cores, 2.4GHz, 1600MHz bus and 16MB L3 cache). Our code was compiled using the g++ 4.8 
compiler with the Cilk Plus extensions for nested parallelism. The only compilation flag we used was the -02 opti- 
mization flag. In all our experiments we use keys of the double data type. The size of each node is about 40 bytes, 
including the two child pointers, the key, the balance information, the size of the subtree, and a reference count. We 
generate multiple sets varying in size from 10 4 to 10 8 . Depending on the experiment the keys are drawn either from 
an uniform or a Gaussian distribution. We use /i and a to denote the mean and the standard deviation in Gaussian 
distribution. Throughout this section n and m represent the two input sizes for functions with two input sets (n > m). 

We test our algorithm by comparing it to other available implementations. This includes the sequential version 
of the set functions defined in the C++ Standard Template Library (STL) [20l and STL’s std: : set (implemented 
by RB tree). The STL supports the set operations set_union, set_intersection, and set.dif ference on 
any container class. Using an std : : vector these algorithms achieve a runtime of 0(m + n). Since the STL does 
not offer any parallel version of these functions we could only use it for sequential experiments. To see how well 
our algorithm performs in a parallel setting, we compare it to parallel WBB-trees fl2l and the MCSTL library 03- 
WBB-trees, as well as the MCSTL, offer support for bulk insertions and deletions. They process the bulk updates 
differently. The MCSTL first splits the main tree among p processors, based on the bulk sequence, and then inserts 
the chunks dynamically into each subtree. The WBB-tree recursively inserts the bulk in parallel into the main tree. To 
deal with heavily skewed sequences they use partial tree reconstruction for fixing imbalances, which takes constant 
amortized time. The WBB-tree has a more cache aware layout, leading to a better cache utilization compared to both 
the MCSTL and our implementation. To make a comparison with these implementations we use their bulk insertions, 
which can also be viewed as an union of two sets. However notice that WBB-trees take the bulk in the form of a sorted 
sequence, which gives them an advantage due to the faster access to the one array than to a tree, and far better cache 
performance (8 keys per cache line as opposed to 1). 


Comparing the balancing schemes and functions To compare the four balancing schemes we choose UNION as 
the representative operation. Other operations would lead to similar results since all operations except JOIN are generic 
across the trees. We compare the schemes across different thread counts and different sizes. 

Figure[TT](a) shows the runtime of Union for various tree sizes and all four balancing schemes across 64 cores. 
The times are very similar among the balancing schemes — they differ by no more than 10%. 

Figure [IT] (b) shows the speedup curves for Union on varying core numbers with two inputs of size 10 8 . All 
balancing schemes achieve a speedup of about 45 on 64 cores, and about 30 on 32 cores. The less-than-linear speedup 
beyond 32 cores is not due to lack of parallelism, since when we ran the same experiments on significantly smaller 
input (and hence less parallelism) we get very similar curves (not shown). Instead it seems to be due to saturation of 
the memory bandwidth. 

We use the AVL tree as the representative tree to compare different functions. Figure[TT](c) compares the Union, 
Intersect and Difference functions. The size of the larger tree is fixed (10 8 ), while the size of the smaller tree 
varies from 10 4 to 10 8 . As the plot indicates, the three functions have very similar performance. 

The experiments are a good indication of the performance of different balancing schemes and different functions, 
while controlling other factors. The conclusion is that all schemes perform almost equally on all the set functions. It 
is not surprising that all balancing schemes achieve similar performance because the dominant cost is in cache misses 
along the paths in the tree, and all schemes keep the trees reasonably balanced. The AVL tree is always slightly faster 
than the other trees and this is likely due to the fact that they maintain a slightly stricter balance than the other trees, 
and hence the paths that need to be traversed are slightly shorter. For different set functions the performance is also as 
expected given the similarity of the code. 
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Figure 11: (a) Times for Union as a function of size ( n = 10 s ) for different BBSTs; (b) speed up of Union for 
different BBSTs; (c) times for various operations on AVL trees as a function of size (n = 10 8 ); (d) comparing STLs 
set union with our Union; (e, f, g, h) comparing our Union to other parallel search trees; (e, h) input keys are 
uniformly distributed doubles in the range of [0, 1]; (f, g) inputs keys follow a normal distribution of doubles - the 
mean of the main tree is always fii = 0, while the mean of the bulk is fi -2 = 1. Figure (f) uses a standard deviation of 
a = 0.25, while Figure (g) shows the performance across different standard deviations. 


Given the result that the four balancing schemes do not have a big difference in timing and speedup, nor do the 
three set functions, in the following experiments we use the AVL tree along with UNION to make comparisons with 
other implementations. 

Comparing to sequential implementations The STL supports set_union on any sorted container class, including 
sets based on red-black trees, and sorted vectors (arrays). The STL set_union merges the two sorted containers by 
moving from left to right on the two inputs, comparing the current values, and writing the lesser to the end of the output. 
For two inputs of size n and m,m<n, it takes 0(m+n ) time on std : : vectors, and 0((n+m) log (n+m)) time 
on std : : set (each insertion at the end of the output red-black tree takes 0(\og(n+m)) time). In the case of ordered 
sets we can do better by inserting elements from the smaller set into the larger, leading a time of 0(mlog(n + m). 
This is also what we do in our experiments. For vectors we stick with the available set _union implementation. 

Figure [XT] (d) gives a comparison of times for Union. For equal lengths our implementation is about a factor of 3 
faster than set variant (red-black trees), and about 8 times slower than the vector variant. This is not surprising since 
we are asymptotically faster than their red-black tree implementation, and their array -based implementation just reads 
and writes the values, one by one, from flat arrays, and therefore has much less overhead and much fewer cache misses. 
For taking the union of smaller and larger inputs, our Union is orders of magnitude faster than either STL version. 
This is because their theoretical work bound (0(m + n) and 0(m log(m + n)) is worse than our 0(m log (n/m + 1)), 
which is optimal in comparison model. 

Comparing to parallel implementations on various input distributions We compare our implementations to other 
parallel search trees, such as the WBB-trees, as described in fl2l . and parallel RB trees from the MCSTL fl3l . We 
test the performance on different input distributions. 
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In Figure [TT| (e) we show the result of Union on uniformly distributed doubles in the range of [0,1] across 64 
cores. We set the input size to n = m = 10', 'i. from 4 to 8. The three implementations have similar performance 
when n = m = 10 4 . As the input size increases, MCSTL shows much worse performance than the other two because 
of the lack of parallelism (Figure 1 1 (h) is a good indication), and the WBB-tree implementation is slightly better 
than ours. For the same reason that STL vectors outperform STL sets (implemented with RB trees) and our sequential 
implementation, the WBB-trees take the bulk as a sorted array, which has much less overhead to access and much 
better cache locality. Also their tree layout is more cache efficient and the overall height is lower since they store 
multiple keys per node. 

Figure 11 (f) shows the result of a Gaussian distribution with doubles, also on all 64 cores with set sizes of 
n = m = 10* for i = 4 through 8. The distributions of the two sets have means at 0 and 1 respectively, and both 
having a standard deviation of 0.25, meaning that the data in the two sets have less overlap comparing to a uniform 
distribution (as in (e)). In this case our code achieves better performance than the other two implementations. For our 
algorithms less overlap in the data means more parts of the trees will be untouched, and therefore less nodes will be 
operated on. This in turn leads to less processing time. 

We also do an in-depth study on how the overlap of the data sets affects the performance of each algorithm. We 
generate two sets of size n = m = 10 8 , each from a Gaussian distribution. The distributions of the two sets have 
means at 0 and 1 respectively, and both have an equal standard deviation varying in {1, 1/2, 1/4, 1/8, 1/16}. The 
different standard deviations are to control the overlap of the two sets, and ideally less overlap should simplify the 
problem. Figure[TT](g) shows the result of the three parallel implementations on a Gaussian distribution with different 
standard deviations. From the figure we can see that MCSTL and WBB-tree are not affected by different standard 
deviations, while our join-based union takes advantage of less overlapping and achieves a much better performance 
when a is small. This is not surprising since when the two sets are less overlapped, e.g., totally disjoint, our Union 
will degenerate to a simple Join, which costs only Of log n) work. This behavior is consistent with the “adaptive” 
property (not always the worst-case) in [?]. This indicates that our algorithm is the only one among the three parallel 
implementations that can detect and take advantage of less overlapping in data, hence have a much better performance 
when the two operated sets are less overlapped. 

We also compare the parallelism of these implementations. In Figure [TT] (h) we show their performance across 64 
cores. The inputs are both of size 10 8 , and generated from an uniform distribution of doubles. It is easy to see that 
MCSTL does not achieve good parallelism beyond 16 cores, which explains why the MCSTL always performs the 
worst on 64 cores in all settings. As we mentioned earlier, the WBB-tree are slightly faster than our code, but when 
it comes to all 64 cores, both algorithms have similar performance. This indicates that our algorithm achieves better 
parallelism. 

To conclude, in terms of parallel performance, our code and WBB-trees are always much better than the MCSTL 
because of MCSTL’s lack of parallelism. WBB-trees achieve a slightly better performance than ours on uniformly 
distributed data, but it does not improve when the two sets are less overlapped. Thus our code is much better than 
the other two implementations on less overlapped data, while still achieving a similar performance with the other 
algorithms when the two sets are more intermixed with each other. 
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7 Conclusions 


In this paper, we study ordered sets implemented with balanced binary search trees. We show for the first time that a 
very simple “classroom-ready” set of algorithms due to Adams’ are indeed work optimal when used with four different 
balancing schemes-AVL, RB, WB ttees and treaps — and also highly parallel. The only ttee-specific algorithm that 
is necessary is the Join, and even the Joins are quite simple, as simple as Insert or Delete. It seems it is not 
sufficient to give a time bound to JOIN and base analysis on it. Indeed if this were the case it would have been done 
years ago. Instead our approach defines the notion of a rank (differently for different Uees) and shows invariants on 
the rank. It is important that the cost of JOIN is proportional to the difference in ranks. It is also important that when 
joining two trees the resulting rank is only a constant bigger than the larger rank of the inputs. This insures that when 
joins are used in a recursive tree, as in UNION, the ranks of the results in a pair of recursive calls does not differ much 
on the two sides. This then ensures that the set functions are efficient. 

We also test the performance of our algorithm. Our experiments show that our sequential algorithm is about 3x 
faster for union on two maps of size 10 8 compared to the STL red-black tree implementation. In parallel settings our 
code is much better than the two baseline algorithms (MCSTL and WBB-tree) on less overlapped data, while still 
achieves similar performances with WBB-tree when the two sets are more intermixed. Our code also achieves 45x 
speedup on 64 cores. 
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A Proofs for Some Lemmas 
A.l Proof for Lemma [8] 

Proof. One observation in WB trees and treaps is that all nodes attached to a single layer root form a chain. This is 
true because if two children of one node v are both in layer i, the weight of v is more than 2 I+1 , meaning that v should 
be layer i + 1. 

For a layer root v in a WB tree on layer k, w(v) is at most 2 k+1 . Considering the balance invariant that its child 
has weight at most (1 — a)w(v), the weight of the i-th generation of its descendants is no more than 2 fc+1 (l — af. 
This means that after t* = log 2 generations, the weight should decrease to less than 2 k . Thus d(v) < log 2, 
which is a constant. 

For treaps consider a layer root v on layer k that has weight N £ [ 2 k , 2 k+1 ). The probability that d(v) > 2 is 
equal to the probability that one of its grandchildren has weight at least 2 k . This probability P is: 
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Figure 12: An illustration of two kinds of outcomes of rotation after joining two weight balanced trees. After we 
append the smaller tree to the larger one and rebalance from that point upwards, we reach the case in (0), where u has 
been balanced, and the smaller tree has been part of it. Now we are balancing v, and two options are shown in (1) and 
(2). At least one of the two rotation will rebalance v. 


N a o k 

E / 

(7) 

t=2 fc + l 


2 k+1 . nk 

E '-=r 

(8) 

i=2 k +l 


In 2 

(9) 


We denote 1 — In 2 as p c . Similarly, the probability that <l(v) > 4 should be less than p 2 , and the probability shrink 
geometrically as d{v) increase. Thus the expected value of d(v) is a constant. 

Since treaps come from a random permutation, all s(v) are i.i.d. □ 


A.2 Proof for Lemma ID 

Proof. We are trying to show that for T r =Union(T p , Td ) on AVL, RB or WB trees, if r{T p ) > r(Td) then r(T r ) < 
r(T p )+r(T d ). 

For AVL and RB trees we use induction on r(T p ) + riTf). When r(Td ) + r(T p ) = 1 the conclusion is trivial. If 
r = r(T p ) > r(Td), T p will be split into two subtrees, with rank at most r(T p ) — 1 since we remove the root. Td will be 
split into two trees with height at most r(Td) (Theorem[2jl. Using the inductive hypothesis, the two recursive calls will 
return two trees of height at most r(T p ) — 1 + r(Td). The result of the final JOIN is therefore at most r(T p ) + r(Td). 
For WB trees, \T\ < \T p \ + \T d \ < 2\T p \. Thus r(T) < r(T p ) + 1 < r(T p ) + r(T d ). □ 
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A.3 Proof for Lemma |3] 

Proof. Recall that in a weight balanced tree, for a certain node, neither of its children is (3 times larger than the other 
one, where (3=^ — 1. When a < 1 — we have (3 > 1 + y/2. 

WLOG, we prove the case when \Tl | < \Tr |, where Tl is inserted along the left branch of Tr. Then we rebalance 
the tree from the point of key k and go upwards. As shown in Figure[l2](0), suppose the rebalance has been processed 
to u (then we can use reduction). Thus the subtree rooted at u is balanced, and Tr is part of it. We name the four trees 
from left to right A, B, C and D, and the number of nodes in them a, b. c and cl. From the balance condition we know 
that A is balanced with B + C, and B is balanced to C, i.e.: 

-p(b + c) <a< (3(b + c) 

^jb <c < Pc 

We claim that at least one of the two operations will rebalanced the tree rooted at v in Figure [T2|(0): 

Op. (1). Single rotation: right rotation at u and v (as shown in Figure [T2|(l)); 

Op. (2). double rotation: Left rotation followed by a right rotation (as shown in Figure [l2|(2)). 

Also, notice that the inbalance is caused by the insertion of a subtree at the leftmost branch. Suppose the size 
of the smaller tree is x, and the size of the original left child of v is y. Note that in the process of JOIN, Tl is not 
concatenated with v. Instead, it goes down to deeper nodes. Also, note that the original subtree of size y is weight 
balanced with D. This means we have: 


GO) 

( 11 ) 


x < d + y ) 

^d<y < pd 
x+y=a+b+c 

From the above three inequalities we get x < jjd + d, thus: 


a + 


b + c= x + y< (1 + /3 + —)d 


Since a unbalance occurs, we have: 

a + b + c > pd 

We discuss the following 3 cases: 

Case 1. B + C is weight balanced with D, i.e., 

~p(b + c) < d< P(b + c ) 


( 12 ) 

(13) 


(14) 


In this case, we apply a right rotate. The new tree rooted at u is now balanced. A is naturally balanced. 

Then we discuss in two cases: 

Case 1.1. (3a > b + c + d. 

Notice that 6 + c > 4 a, meaning that b+ c+d > j^a. Then in this case, A is balanced to B + C+D, 
B + C is balanced to D. Thus just one right rotation will rebalance the tree rooted at u (Figure [l2| 

(D). 
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Case 1.2. [3a < b + c + d. 

In this case, we claim that a double rotation as shown in Figure[l2](2) will rebalance the tree. Now 
we need to prove the balance of all the subtree pairs: A with B, C with D , and A + B with C + 1). 
First notice that when /3a < b + c + d, from we can get: 

f3d<a + b + c< — ( b + c+d) + b + c 

=K0 — -p)d < (— + 1)(6 + c) 

=>(f3 — l)d < b + c (15) 


Considering ( 14 1 , we have (/? — 1 )d < b + c < /3d. Notice b and c satisfy ( 1 1 1 , we have: 

b 


0 — 1 


> 3+T (i, + 0)> ^TT‘ 


c> fltT ( ' , + <;)> fTT <i 


(16) 

(17) 


Also note that when /? > 1 + y/2 « 2.414, we have 


0 + 1 
0-1 


<0 


We discuss the following three conditions of subtrees’ balance: 

I. Prove A is weight balanced to B. 

i. Prove b < [3a. 

Since [3a < b + c (applying , we have b < [3a. 

ii. Prove a < /3b. 

In the case when [3a < b + c + d we have: 

a < ^ (6 + c + d) (applying (fll|) , < |l6| ) ) 

_ 0 + 1 r 


< 06 


(18) 


II. Prove C is weight balanced to D. 

i. Prove c < (3d. 

Since b + c < 0d (applying fT4|)< we have c < 0d 

ii. Prove d < /3c. 

From ([l7]>, we have 


d < 


0 + 1 

0-1 


c < 0c 


III. Prove A + B is weight balanced to C + 13. 
i. Prove a + b < /3(c + d). 
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From CD.GD and ( p~8j ) we have: 

1 


d< 


0-1 
fi + 1 


(b + c) < 


- 1 


(/3c + c) 


0-1 
>4d < c 


c < / 3c 


=Ki + < (i + 0) c 

=>(1 + ^ + P)d < 0(c + d) + c (applying ( [l2| ) 

=^0 + 6 + c< (1 + — + 0)d < /3(c + d) + c 
=>a + b < fi{c + d) 

ii. Prove c + d < fi(a + b). 

When /3 > 2, we have < 0- Thus applying ( 15 i and ( 10 1 we have: 

d < -^-j-(d + c) < < 0 a 


Also we have c < fib (applying (111). Thus c + d < 0(a + b). 


Case 2. B + C is too light that cannot be balanced with D, i.e., 

0(6 + c) < d 


(19) 


In this case, we have a < 0(6 + c) < d (applying ( 10 1 and ( 19 1), which means that a + b + c<d+jjd< fid 

when fi > 1+ 2 V ^ ss 1.618. This contradicts with the condition that A + 33 + C is too heavy to D (a + b + c > 
fid). Thus this case is impossible. 


Case 3. B + C is too heavy that cannot be balanced with D , i.e., 

b+ c> fid 
=>a > — ( b + c ) > d 


( 20 ) 

( 21 ) 


In this case, we apply the double rotation. 

We need to prove the following balance conditions: 


I. Prove A is weight balanced to B. 


i. 

ii. 


Prove b < /3a. 

Since fia > b + c (applying ( 10 1 ) , we have b < fia. 

Prove a < fib. 

Suppose c = kb, where ^ < k < fi. Since b + c > fid, we have: 


d < 



( 22 ) 
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From the above inequalities, we have: 


\ + b + c = a + b + kb ( applying (fl2|)) 

<(l+/3+-)d (apply ing©) 


P 

<(! +P + ^)x—b 


P 


1 


a < 


-1 )(1 + k)b 


^(i+w 


< 


(P+ l ) 2 
P 2 


When/? > I x ( ^y+ 47 ) ~ 1/3 + ( ^y+ 47 ) : 1 / 3 + 1 « 2.1479, we have < /?. Hence a < /3&. 

H. Prove C is weight balanced to D. 

i. Prove d < /3c. 

When /3 > 1+ 2 v/ ^ 1.618, we have p > 1 + j. Assume to the contrary c < jjd, we have 

b < Pc < d. Thus: 

b + c < (1 + —)d < pd 

, which contradicts with (|20| that B + C is too heavy to be balanced with D. 

ii. Prove c < (3d. 

Plug (21 1 in ( 12 1 and get b + c < (P + jj)d. Recall that /3 > 1, we have: 

— c + c < b + c < (P + — )d 

p 2 + 1 , Q , 

=>C < p + 1 d < pd 

III. Prove A + B is weight balanced with C + D. 

i. Prove c + d < (3(a + b). 

From ( 1 1 1 we have c < Pb, also d < a < Pa (applying (21 1 ), thus c + d < P(a + b). 

ii. Prove a + b < (3{c + d). 

Recall when P > 1+ ^ & 1.618, we have /? > 1 + j. Applying (20 1 and jllj we have: 

d< ^-(b + c) < c + ^-c < Pc 

P p 

^ ~p d< ° 

=t>(l + —)d < (1 + P)c 

=t>(l + + P)d < P(c + d) + c (applying @) 

=$*a ~F 6 4- c < (1 4“ — T fi)d < p{c 4~ (?) 4~ c 
=>a + b < P(c + d) 
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Taking all the three conclusions into consideration, after either a single rotation or a double rotation, the new subtree 
will be rebalanced. 

Then by induction we can prove Lemma[3] □ 
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